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POLARIZATION 


QUARTZ 


1 < a <30 
n~a -4 



0 20 40 60 80 100 120 140 160 180 

SCATTERING ANGLE (DEG) 

Figure A-15. Comparison of the degree of polarization for quartz 
obtained from the trapezoidal rule and integration 
increments of A a = 0.1, 0.3, and 0.5. 
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Figure A-14. Total intensity curves for quartz as obtained 
from the trapezoidal rule and integration increments 
of A o; = 0.1, 0.3, and 0.5. 





POLARIZATION 
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INDEX OF REFRACTION m « 1.33-.00i 
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SIZE PARAMETER RANGE a mjn - 1 TO a mul - 30 
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5 POINT GAUSS-LEGENDRE 



SCATTERING ANGLE (DEG) 

Figure A-13. A comparison of the degree of polarization 
from three Gauss -Legendre integration schemes . 
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Figure A-ll. Total scattering intensity function from the trapezoidal rule 
for a mixture of half quartz and half ice (A a = 0. 1) . 
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Figure A-10. Degree of polarization for 50 percent quartz, 50 percent ice 
mixture using the 96 -point Gaussian quadrature 
(a min = 1.0, a max = 20. 0). 
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Figure A-9. Total scattering intensity function for 50 percent quartz, 
50 percent ice using the 96-point Guassian quadrature 
{a. min = 1.0, a max = 20.0). 
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ZINC 

M = 1.41 - 4, 10i 
1 £ a s 30 



scattering angle 'degrees! 

“96- POINT GAUSS -LEGENDRE INTEGRATION 


Figure A-8. Degree of polarization for zinc using 
the 96 -point Gauss -Legendre integration . 
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ZINC 

M = 1.41 - 4. lOi 

1 =S a 5 30 
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3CA TTEfil NC ANCLE <CECrEES> 

5-POINT GAUSS- LEGENDRE INTEGRATION 


Figure A -7. Degree of polarization for zinc from 
the 5-point Gauss-Legendre integration . 
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ZINC 

M = 1.41 - 4. lOi 
1 £ a ^ 30 



scattering ANGLE (CEGREES) 

96-POINT GAUSS -LEGENDRE INTEGRATION 

Figure A-6. Total scattering intensity function for zinc using 
the 96 -point Gauss -Legendre integration . 
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ZINC 

M = 1.41 - 4. lOi 
1 < a s 30 



SCATTERING AnGcE tCEGREES) 

5- POINT GAUSS- LEGENDRE INTEGRATION 


Figure A-5. Total scattering intensity function for zinc using 
the 5-point Gauss-Legendre integration . 
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WATER (ICE) 

M = 1.33 - O.OOi 
1 < a < 30 



0 20 4 0 6 0 60 100 IPO M'l * ' tf-0 1 60 TOO 

*'< a* V t'l(. AU.t \ ‘ t t 


96- POINT GAUSS- LEGENDRE INTEGRATION 


Figure A-4. Degree of polarization for water using 
the 96 -point Gauss-Legendre integration . 
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WATER (ICE) 

M = 1.33 - O.OOi 
1 < a s 30 



SCATTERING ANGLE (CE&rEES) 

5- POINT GAUSS- LEGENDRE INTEGRATION 

Figure A-3. Degree of polarization for water 
from the 5-point Gauss-Legendre integration . 
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Figure A-2. Total scattering intensity function for water using 
the 96 -point Gauss -Legendre integration . 
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WATER (ICE) 

M = 1.33 - O.OOi 
1 < a s 30 



SCATTERING ANCLE (DECREES' 

5- POINT GAUSS- LEGENDRE INTEGRATION 


Figure A-i. Total scattering intensity function for water from 
the 5-point Gauss-Legendre integration. 
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TABLE A-l. COMPARISON OF 96-POINT GAUSS-LEGENDRE 
AND TRAPEZOIDAL QUADRATURE FOR MIXTURE 
OF QUARTZ AND ICE (l < a < 20) 


Scattering 

Angle 

Gauss-Legendre 

Trapezoidal Quadrature 

Intensity 

Polarization 

Intensity 

Polarization 

KM 

27.517848 

0.000000 

27.467983 

0.000000 


21.930220 

-0.001570 

21.912843 

-0.001516 

10.00 

13.933973 

-0.003175 

13; 936104 

-0.003217 

15. 00 

9.627869 

-0.002480 

9.627133 

-0.002501 

20. 00 

6.968873 

. -0.000947 

6.968865 

-0.000810 

25.00 

5. 170425 

0.003459 

5.169204 

0.003331 

30.00 

3.881988 

0.009044 

3.881645 

0.009103 

35. 00 

2.951191 

0.017750 

2.950012 

0.017596 

40. 00 

2.261300 

0.028241 

2.260898 

0. 028135 

45.00 

1.749777 

0.042031 

1.749769 

0.041972 

50.00 

1.365778 

0.058184 

1.365372 

0.058073 

55. 00 

1.074961 

0.077380 

1.075145 

0.076950 

60.00 

0.855332 

0.098211 

0.854572 

0.097857 

65.00 

0.685607 

0.120485 

0.685912 

0.120480 

70.00 

0.556637 

0.143266 

0.555991 

0.142673 

75.00 

0.456271 

0.164757 

0.456058 

0.164959 

80.00 

0.378865 

0.184513 

0.378550 

0. 184058 

85.00 

0.319032 

0.199275 

0.318594 

0.200081 

90.00 

0.272443 

0.210949 

0.272230 

0.210609 

95.00 

0.236955 

0.213049 

0.236601 

0.214309 

100.00 

0.209769 

0.210831 

0.209368 

0.210396 

105.00 

0. 189428 

0. 195691 

0.189207 

0.196939 

110.00 

0.174989 

0.174005 

0.174603 

0.173860 

115.00 

0.165090 

0. 140950 

0.164921 

0.141842 

120.00 

0. 159906 

0.100688 

0.159734 

0.100845 

125.00 

0. 158751 

0.054130 

0.158466 

: 0. 054847 

130.00 

0. 161664 

0.005797 

0.161531 

0.006187 

135.00 

0. 169249 

-0. 042521 

0.168930 

-0.042707 

140.00 

0.180715 

-0.087857 

0.180417 

-0.087401 

145.00 

0. 197115 

-0.128472 

0. 197125 

-0.127842 

150.00 

0.217899 

-0.165958 

0.217781 

-0.166240 

155.00 

0.242010 

-0.198321 

0.241537 

-0.199711 

160.00 

0.267969 

-0.227409 

0.268076 

-0.227614 

165.00 

0.298947 

-0.207388 

0.298369 

-0.206938 

170.00 

0.357341 

-0. 167437 

0.355763 

-0.173064 

175.00 

0.422300 

-0.079678 

0.420986 

-0.079102 

180.00 

0.515393 

0.000000 

0.509820 

0.000000 
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the 5- and 96-point integration. All four figures are for water m = 1.33 
-O.OOi and a particle size range of 1 ^ a ^ 30. The same effect can be 
seen in Figures A-5 through A-8, which are for the same integrations except 
that a metal (zinc, m = 1.41 - 4.10i) has been used as the material. The 
smoothing effect and the dominance of the small particles is quite noticeable, 
particularly in the degree of polarization graphs. Figures A-9 through A-12 
compare the results of the integration of the 96 -point Gauss-Legendre and the 
trapezoidal quadrature with A a = 0.1. A mixture of 50 percent quartz and 
50 percent ice particles was used in these graphs. Since the graphs are not 
too accurately plotted, Table A-l, which gives the values obtained by each 
quadrature, is included so that a better comparison can be made. 

Figure A -13 plots the degree of polarization for water that results 
from each Gauss-Legendre quadrature. The convergence toward one smooth 
curve is obvious. Figure A-14 shows the total intensity obtained using the 
trapezoidal rule and various A a' s. It should be noted that this curve is 
insensitive to the integration increment (up to A a = 0.5 ). Figure A-15 
plots the degree of polarization that results from the trapezoidal quadrature 
for various A a' s. These results agree entirely with those of Dave [7] , who 
used a distribution of the form 


n(aj = a a e 

He concludes that A a = 0. 1 is the largest size interval that can be used to 
obtain reliable results. With size intervals of A a = 0.5 and 1.0, ripples 
appear in some angular regions. For A a = 2.0 pseudo-features appear in 
the curves. (Compare Figs. A-3 and A-7 with Figs. A-4 and A-8, respec- 
tively. Also see Fig. A-13. ) Such pseudo-features would tend to invalidate 
any results from model calculations in which the scattering functions obtained 
from such a coarse integration were used. 


132 



DPTT'W)TMr! T>bm? PT AW VnT PTTA/rrr 


APPENDIX 

RESULTS OF THE STUDY 
OF INTEGRATION SCHEMES AND INCREMENTS 

As noted previously, two integration schemes were used, a Gauss- 
Legendre and a trapezoidal quadrature. It was found that large integration 
increments give higher values of intensity than resulted from the small 
increments. The polarization curves also differed considerably. For the 
particle size range of 1 ^ s 30, a 16-point Gauss-Legendre scheme was 
originally used. However, this was later modified to a 96 -point Gauss- 
Legendre scheme and then to one using either the 16- or 96 -point but restrict- 
ing the integration limits so that oi -a . = 1.6 or 9.6, respectively, 

max min 

and changing the upper and lower limits, i. e. , for 1 ^ o’ ^ 5 for the 
16 -point scheme 


1.6 3.2 

/ i<». a, 9)n(a)da + f i(m, a, 0)n(o')dQ! 

1 1.6 


4.8 

+ f i(m, a, 0)n(o;)da 
3.2 


5.0 

+ J i(m, a, 0)n(a)da’ 
4.8 


For the computer program, to reduce the number of integrations needed, any 

excess, such as occurs in the fourth term in the above equation, was placed 

in the first integral. In no case was the integration interval allowed to exceed 

A a = 0.2. It should be noted that the Gauss-Legendre schemes do not work 

well when a is larger than the number of zeros of the Legendre polyno- 
max 

mial (cos d). In this situation the integration is weighted disproportion- 
ately toward the larger particles, negating the smoothing effect the integration 
should have. 

Figures A-l and A-2 give the results for total radiance for water 
(m = 1.33 - O.OOi) using a 5- and 96 -point Gauss-Legendre integration. 
Figures A-3 and A-4 compare the degree of polarizations that results from 
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ZINC 

M = 1.41 - 4.101 

1 s a s 25 



SCA TTZP! N6 SKSU (OEfcjKsi 


Figure 81„ Polarization versus scattering angle 
for zinc with 1 s a =£ 25. 
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ZINC 

M = 1.41 - 4.101 



Figure 80. Polarization versus scattering angle 
for zinc with 1 — a ^ 20. 



ZINC 

M = 1.41 - 4. lOi 

1 < a £ 15 



Figure 79. Polarization versus scattering angle 
for zinc with 1 — a — 15. 


126 




Figure 78. Polarization versus scattering angle 
for zinc with 1 ^ a. — 10. 
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ZINC 

M = 1.41 - 4. 10i 



SCATTERING ANGLE (0EGREE3) 


Figure 77. Polarization versus scattering angle 
for zinc with 1 ^ o; < 5. 
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IRON 

M = 1.27 - 1.371 



Figure 76. Polarization versus scattering angle 
for iron with 1 ^ a — 25. 
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IRON 

M = 1. 27 - 1. 37i 
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Figure 75. Polarization versus scattering angle 
for iron with 1 — a ^ 20. 
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IRON 

M = 1.27 - 1. 37i 



Figure 74. Polarization versus scattering angle 
for iron with 1 — oi — 15. 
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IRON 

M = 1. 27 - 1. 37i 
1 £ a £ 10 



scattering anh.e<oecRee3> 


Figure 73. Polarization versus scattering angle 
for iron with 1 ^ a s 10. 
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SCATTERING ANGLE (DEGREES) 


Figure 72. Polarization versus scattering angle 
for iron with 1 < a < 5. 
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Figure 71. Polarization versus scattering angle 
for graphite with 1 ^ a ^ 25. 
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GRAPHITE 
M = 1.59 - 0. 66i 
1 < a < 20 



Figure 70. Polarization versus scattering angle 
for graphite with 1 < a ^ 20. 
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GRAPHITE 
M = 1. 59 - 0. 66i 
1 £ a s 15 



Figure 69. Polarization versus scattering angle 
for graphite with 1 ^ a £ 15. 
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GRAPHITE \ 

M = 1.59 - 0.661 

1 s Q < 10 ! 
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Figure 68. Polarization versus scattering angle 
for graphite with 1 ^ a ^ 10. 
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SCATTERING ANGLE (DEGREES* 


Figure 67. Polarization versus scattering angle 
for graphite with 1 < a ^ 5. 
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WATER (ICE H) 
M = 1.33 - 0. lOi 



Figure 66. Polarization versus scattering angle 
for water (ice II) with 1 < a < 25. 
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WATER (ICE H) 

M = 1.33 - 0. 10 i 
1 s a < 20 



«o «o to 100 ito uo i«o r 5 o too 
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Figure 65. Polarization versus scattering angle 
for water ( ice II) with 1 < a < 20. 
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WATER (ICE R) 

M = 1.33 - 0. lOi 



Figure 64. Polarization versus scattering angle 
for water (ice II) with 1 < a £ 15. 
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Figure 63. Polarization versus scattering angle 
for water (ice II) with 1 < q? < 10. 
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WATER (ICE R) 

M = 1.33 - 0. 10 i 
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Figure 62. Polarization versus scattering angle 
for water ( ice II) with 1 < a ^5. 


109 


WATER (ICE I) 

M = 1.33 - 0.05i 
1 s a S 25 
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Figure 61. Polarization versus scattering angle 
for water ( ice I) with 1 < a < 25. 
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WATER (ICE I) 

M = 1.33 - 0.05i 



Figure 60. Polarization versus scattering angle 
for water ( ice I) with 1 < - a ^ 20. 
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i | WATEH (ICE I) 

! M = 1.33 - 0.05i 
! 1 s a s 15 



Figure 59. Polarization versus scattering angle 
for water (ice I) with 1 < a < 15. 
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Figure 58. Polarization versus scattering angle 
for water ( ice I) with 1 — oi — 10. 
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WATER (ICE I) 

M = 1.33 - 0.05i 



Figure 57. Polarization versus scattering angle 
for water ( ice I) with 1 < a < 5. 
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TITANITE 
, M = 2. 00 - 0. OOi 
! 1 £ Oi £ 25 



V. Of; 

Figure 56. Polarization versus scattering angle 
for titanite with 1 — a £ 25. 
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TITANITK 
M = 2. 00 - 0. OOi 
1 £ a < 20 



Figure 55. Polarization versus scattering angle 
for titanite with 1 — a — 20. 
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TITAN1TE 

M = 2.00 - 0.001 
1 s a < 15 



Figure 54. Polarization versus scattering angle 
for titanite with 1 £ a ^ 15. 
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TITANITE 
M = 2.00 - 0. OOi 
Isas 10 




TITANITE 
M = 2. 00 - 0. 001 
i £ O' £ 5 



S ' A T T r M no a HGI E ir.c C C CE S » 

Figure 52. Polarization versus scattering angle 
for titanite with 1 ^ a < 5. 
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Figure 51. Polarization versus scattering angle 
for water (ice) with 1 — a — 25. 
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WATER (ICE) 

M = 1.33 - 0. OOi 
1 s a < 20 
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Figure 50. Polarization versus scattering angle 
for water (ice) with 1 < « < 20. 
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> WATER (ICE) 


M = 1.33 - O.OOi 
1 =s a 15 



Figure 49. Polarization versus scattering angle 
for water (ice) with 1 ^ a ^ 15. 
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WATER (ICE) 

M = 1.33 - O.OOi 
1 £ a s 5 



Figure 47. Polarization versus scattering angle 
for water (ice) with l < a s 5 , 




Figure 46. Polarization versus scattering angle 
for quartz with 1 - a — 25. 
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QUARTZ 

I M = 1.55 - O.OOi 

1 < a ST 20 
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Figure 45. Polarization versus scattering angle 
for quartz with 1 - a - 20. 
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scmtcrinc ancle (OeorEEsi 


Figure 44. Polarization versus scattering angle 
for quartz with 1 — a ^ 15. 


91 


POLARIZATION 



90 




QUARTZ 

M = 1. 55 - O.OOi 



Figure 42. Polarization versus scattering angle 
for quartz with 1 ^ a -- 5. 
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ZINC 

M = 1.41 - 4. lOi 
1 £ a < 25 



Figure 41. Intensity versus scattering angle 
for zinc with 1 — a s 25. 
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ZINC 

M = 1.41 - 4. lOi 
1 S a < 15 



Figure 39. Intensity versus 
for zinc with 1 — a 


scattering angle 
£ 15. 



ZINC 

M = 1.41 - 4. lOi 
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Figure 38. Intensity versus scattering angle 
for zinc with 1 s a < io, 
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ZINC 

M = 1.41 - 4. lOi 
1 ^ a < 5 
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SCATTERING ANGLE (DEGREE ST 


Figure 37. Intensity versus scattering angle 
for zinc with 1 < a ^ 5. 



IRON 

M = 1.27 - 1. 37i 
1 s a s 25 



Figure 36, Intensity versus scattering angle 
for iron with 1 - a - 25, 
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Figure 35. Intensity versus scattering angle 
for iron with l < a < 20. 
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IRON 

M = 1.27 - 1. 37i 
1 s a s 15 



Figure 34. Intensity versus scattering angle 
for iron with 1 ^ a ^ 15. 
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IRON 

M = 1.27 - 1. 37i 
1 < a s io 
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Figure 33. Intensity versus scattering angle 
for iron with 1 ^ 0! £ 10. 
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IRON 

M = 1.27 - 1. 37i 

1 s Q! < 5 



SCATTERING ANGLE (DEGREES) 


Figure 32. Intensity versus scattering angle 
for iron with 1 ^ a < 5. 



GRAPHITE 
M = 1. 59 - 0. 66i 
1 s a s 25 



scttcimi. ««,le (OE(.aee>i 


Figure 31. Intensity versus scattering angle 
for graphite with 1 ^ a ^ 25. 
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GRAPHITE 
M = 1.59 - 0. 66i 
1 S a < 20 



Figure 30. Intensity versus scattering angle 
for graphite with 1 £ a <20. 
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GRAPHITE 
M = 1. 59 - 0. 66i 
1 < a ^ 15 



Figure 29. Intensity versus scattering angle 
for graphite with 1 < a ^ 15. 



GRAPHITE 
M = 1.59 - 0.661 



Figure 28. Intensity versus scattering angle 
for graphite with 1 :£ a ^ 10; 
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Figure 27. Intensity versus scattering angle 
for graphite with 1 ^ a. £ 5. 





WATER (ICE H) 

M = 1.33 - O.lOi 
1 £ a 25 



scATU«iNo *Ntri_e 


Figure 26. Intensity 
for water ( ice II) 


versus scattering angle 
with 1 < a < 25. 
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SCATTCRIHO 


Figure 25. Intensity ver; 
for water ( ice II) wil 


water (ice n) 

M = 1.33 - 0.11 
1 ^ a =S 20 


3 



scattering angle 
< o; < 20. 






WATER (ICE H) 
M = 1. 33 - O.lOi 
1 s a s 15 




Figure 24. Intensity versus scattering angle 
for water (ice II) with 1 < a <15. 
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WATER (ICE E) 

M = 1.33 - O.lOi 
1 s a < 10 
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Figure 23. Intensity versus scattering angle 
for water (ice II) with 1 < a < 10. 





WATER (ICE D) 
M = 1.33 - 0. lOi 
1 s a < 5 



Figure 22. Intensity versus scattering angle 
for water (ice II) with 1 < a — 5. 
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WATER (ICE p 
M = 1.33 - 0.051 





Figure 21. Intensity versus scattering angle 
for water ( ice I) with 1 < a <25. 
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WATER (ICE I) 

M = 1.33 - 9.051 
1 < a < 20 



Figure 20. Intensity versus scattering angle 
for water ( ice I) with 1 < « < 20. 
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Figure 19. Intensity versus scattering angle 
for water ( ice I) with 1 < a < 15. 
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WATER (ICE I) 

M = 1.33 - 0.05i 
1 £ a =S 10 



Figure 18. Intensity versus scattering angle 
for water (ice I) with 1 < q < 10. 
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WATER (ICE I) 

M = 1.33 - 0. 05i 
1 < tt S 5 



Figure 17. Intensity versus scattering angle 
for water ( ice I) with 1 < a < 5. 
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TITANITE 

M = 2.00 - O.OOi - 
1 s a s 25 





Figure 16. Intensity versus scattering angle 
for titanite with 1 — ca ^25. 



TITANITE 
M = 2.00 - O.OOi 
1 £ a s 20 



Figure 15. Intensity versus scattering angle 
for titanite with 1 ^ a — 20. 
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TITANITE 
M = 2.00 - O.OOi 
1 ^ a s 15 



Figure 14. Intensity versus scattering angle 
for titanite with 1 ^ a ^ 15, 
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TITANITE 
M = 2. 00 - O.OOi 
1 < a ^ 10 



Figure 13. Intensity versus scattering angle 
for titanite with 1 < a. s 10* 
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TITANITE 
M = 2.00 - 0.001 
1 £ a s 5 



Figure 12. 
for 


Intensity versus scattering angle 
titanite with 1 ^ a ~ 5. 
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WATER (ICE) 

M = 1.33 - O.OOi 
1 i< a ^ 25 
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Figure 11. Intensity versus scattering angle 
for water (ice) with 1 ^ a ^ 25. 
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WATER (ICE) 

M = 1.33 - O.OOi 
1 < u< 20 
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Figure 10. Intensity versus scattering angle 
for water (ice) with 1 — a — 20. 
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WATER (ICE) 

M = 1.33 - O.OOi 
1 — tt — 15 



Figure 9. Intensity versus scattering angle 
for water (ice) with 1 s a < 15 . 





WATER (ICE) 

M = 1.33 - O.OOi 
1 < a ^ 5 



Figure 7. Intensity versus scattering angle 
for water (ice) with 1 — a — 5. 



QUARTZ 

m = 1. 55 - 0. OOi 
1 s a < 25 



Figure 6. Intensity versus scattering angle 
for quartz with 1 — a — 25. 
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QUARTZ 

m = 1. 55 - O.OOi 
1 s a < 20 
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Figure 5. Intensity versus scattering angle 
for quartz with 1 — a ^ 20. 
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QUARTZ 

M = 1.55 - O.OOi 
1 £ a <15 



Figure 4. Intensity versus scattering angle 
for quartz with I s a s 15. ' 
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INTENSITY 1 + INTENSITY 2 



0 20 40 60 30 100 120 140 160 180 

SCATTERING ANGLE (DEGREES) 

Figure 3. Intensity versus scattering angle 
for quartz with 1 ^ a — 10. 


50 



QUARTZ 

M = 1. 55 - O.OOi 

1 s a < 5 



SCaMEM*** AN&'-ClQCC-^CE 3 i 


Figure 2. Intensity versus scattering angle 
for quartz with 1 — a < 5. 
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TABLE 40. INTENSITY AND POLARIZATION FOR ZINC (l 
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TABLE 39. INTENSITY AND POLARIZATION FOR ZINC (l 
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TABLE 38. INTENSITY AND POLARIZATION FOR ZINC (l 
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TABLE 37. INTENSITY AND POLARIZATION FOR ZINC (l 
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TABLE 36. INTENSITY AND POLARIZATION FOR ZINC (l 
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INTENSITY AND POLARIZATION FOR IRON (l 
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TABLE 34. INTENSITY AND POLARIZATION FOR IRON (l 
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TABLE 33. INTENSITY AND POLARIZATION FOR IRON (l 
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TABLE 32. INTENSITY AND POLARIZATION FOR IRON (l 
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TABLE 31. INTENSITY AND POLARIZATION FOR IRON (l 


LO 

VI 

3 

VI 


o 

o 

CD 


__ . . . _ . _nNNinffl'0-fflOi>-ocato oxccouminrixfNXxo 

O © --3X X 9 O O O O J X OIVN f .)fM © > «s IN — 

r T 3 - N O 3 N ^ ^ D O O N O .1 ^ -n T 3 0 -I> O Jfl 

ON«. , )nN0--OTT-0‘NS-®TNl/niOai>} > TCD'0-0»(*l-Or-N 
O©<Na0«N©xxx - r\ N N(N -0 — x o x © o 0 8 <0 T N -* Q 
r T J) J) « © -0-0 -O ^3 ifl T r ^^^----0 3 3 0 :0 O 


o 

© 

IN © 

IN 

X 

T 

iD 

o 

© 

— 43 

O' 

X 

9 

O 

•3 

O 

x *> 

f3 

X 

<N 

o 

o 

x 

— '.*3 

o 

X 

O' 

mm 

3 

O 

O rs 

IN 

3) 

Ul 

O 

3 

3 

3 '.3 


— 

IN 

O 


X 

< 

CL 


N lll»N —» CD X X l/l ©XX -0 0-0 3 0 tj <© O X J1 X CO X O T J) N OX 3 N -O N N 

o o > J - < ifl n ■o^iP'Cjo-ajNn^a o ^ <o co co in o 3 in •« a o <i ui a ^ 

i o ~ x x co cmco o oo o- -• -h x j-N-o^raji-N 3 o 3 — r — co « n r ®n 

X 2 ^ lA N B N J N r © X O :N > 3-«'Ofl JNili/IN !) r O' ^ N N IN- IN o 

X 3® f 3 O' » T N 3N O' i/t ifl N 3-OfNX©XO — 3>®®N-0'0'0-0'OB<0® 
O IN © X — rsl r J®ilQ®K4)ininT T'n(nnt3<n«*)C^NNNNNMNNr^NN 

© © . Ji T X O IN N — — — 


a 

: Z 

X 

iD 

O x 

o 

o 

x 

3 

X 

o 

3 0 2 

© 

CO 

o 

IN 

9 

rnm 

IN O X 

a> 

a> 

rs. 

in 

© 

— in 

•3 

X 

in 

<N 

IN 

X — 

o 

a* 

0 

x 

CO JO 

o 

X 

IN 


x 

o 

— 3«J1 

X 

■X 

9 

-O 

CO 

-o 

© x 

© 

X 


9 

rs 

© 

3 O 

(O 

•m 

co 

3 <N 

rs in 

a - 

. , ■ . 

O 

<N 

x m 

<N 

in 

O 

mm 

CO 

X 

x in in 

o 

x 

n 

in 

O 

a 

o m 

>o 

X 

IN 

© 

o 

X 

© rs 

a 

rs 

m 

m 

•M 

x m 

x 

z . 

X 

O* 

o -0 

O 

x 

in 

in 

rm 

3 

J m - 


X 

x 

mm 

IN 

x 

X IN 

•m 

3 

tO 

X 

O' 

X 

O IN 

rs 

■mm 

© 

mm 

© 

X IN 

m 

— • 

O' 

iD 

— O' 

o 

x 

IN 

>0 

X 

Q 

-s x ® 

IN 


© 

© 

© 

-0 

X JO 

O 

O 

3 

— • 

-• 

IN 

IN tO 

o 

X 

X 

in 

in 

in © 

• 


o 

— 

m in 

in 

© 

o 

o 

o 

X 

in m in 

IN 

— 

-• 


— 

— 

*• — 

— 

IN 

IN 

IN 

IN 

«N 

IN IN 

IN 

IN 

IN 

IN 

IN IN IN 


in 

■»< 


mj 

< 


3-,. • 

■ - X 

2 

rs n 

o o 

_ 

cO 

3 

rs 

rs 


9 

8 

rs 

O 

_ m 

X 

n 

© 

X 

n 

•N 

.1 

8 

«• 

m 

N 

> 

X 

X 

ji 

O 

o 

X 

3 O 

3 

3 

i!l 

UJ 

-3 co 

“• *1 

X 

r» 

IN 

N 

— « 

O 

IN 

> 

rs 

-m 

3 rs 

© 

X 

"O 

X 

3 


IN 

X 

X 

X 

X 

IN 

IN 

— 

N 

*• 

3 

IN X O 

© 

3 

. ! *J 

r- 

O X 

in in 

X 

X 

rs 

•0 

X 

X 

X 

IN 

IN 

8 

o co 

CO 

X 

X 

9 

CO 

»n 

X 

IN 

3 

© 

© 

© 

<N 


3 

X 

3 

X 

X 1 X 

rs 

.. 

z 

. x © 

CO .— 

IN 

© 

IN 

X 

X 

in 

rs 

rs 

rs 

© 

X !» 

CO 

O 

X 

CO 

3 x 

o 

O 

X 

X 

in 

X 

© 

in 

© 

CO 

O 

to 

in x 

O 

IN 

2 

—« 

x n 

© X 

•M 

«• 

A. 

~m 

o 

mm 

"1 

m» 

to 

to 

n in 

n 

co 

«■ 

J1 

— • 

0 

o 

3 

X 

in 

o 

mm 

'D 

X 

«© 

X 

X 

© 

© •© 

© 

• 

‘ *- 


O IN 

co in 

© 

X 

IN 

rs 

IN 

X 

© 

X 

IN 

x to 

rs 

© 

© 

in 

in 

X 

X 

X 

o 

o 

O 

o 

o 

IN 

IN 

IN 

IN 


IN IN 

IN 


2 

— 

© © 

in in' 

'X 

o 

>0 

IN 

IN 

J* 

L 

— 

•1 

— 























UJ 

O' 

< n 
**0 
? V/ 
ui z y 

WOw, 

Ul qc yi 
* — '**■ 
CL 

in 


3 3 3 3 3 J 3 I 3 3 3 J 3 1 3 3 3 3 3 3 3 3 3 3 3 3 3 J 3 3 3 3 3 1 3 )3 
O OO'J O JO O 3 3 3 3/3 *.3 333333333033330030 D 030 H 


J }j> )in 3 m 3 in -3 in J J> j A 3 Jl 3 J1 3 J1 )JI 3 J> 
-*-<NNn'*l®T l AiA®'0NN«®>O' 3 3 ■< «* N 


.3 ifl -jifl :i Ji 'iifl-aiii *i 

x in irt « •« nn « 


UJ 

X 


38 


263H^U • 26.3*t60 •OUOoQO 



TABLE 30. INTENSITY AND POLARIZATION FOR GRAPHITE (l 
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TABLE 29 . INTENSITY AND POLARIZATION FOR GRAPHITE (l 
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TABLE 28. INTENSITY AND POLARIZATION FOR GRAPHITE ( 
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TABLE 27. INTENSITY AND POLARIZATION FOR GRAPHITE (l 
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TABLE 26. INTENSITY AND POLARIZATION FOR GRAPHITE (l 
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TABLE 24. INTENSITY AND POLARIZATION FOR ICE II ( 1 < a 
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TABLE 23. INTENSITY AND POLARIZATION FOR ICEn(l< a < 15) 
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TABLE 22. INTENSITY AND POLARIZATION FOR ICE II ( 1 < a < 10) 
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TABLE 21. INTENSITY AND POLARIZATION FOR ICE II ( 
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TABLE 20. INTENSITY AND POLARIZATION FOR ICE I ( 1 a ^25) 
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TABLE 19. INTENSITY AND POLARIZATION FOR ICE I ( 1 < a 
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TABLE 18. INTENSITY AND POLARIZATION FOR ICE I ( 1 < a 
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TABLE 17. INTENSITY AND POLARIZATION FOR ICE I ( 1 < 
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TABLE 16. INTENSITY AND POLARIZATION FOR ICE I (1 < 
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TABLE 15. INTENSITY AND POLARIZATION FOP TITANITE (l 
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TABLE 14. INTENSITY AND POLARIZATION FOR TITANITE (l 
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TABLE 13. INTENSITY AND POLARIZATION FOR TITANITE (l 
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TABLE 12. INTENSITY AND POLARIZATION FOR TITANITE (l * a 
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TABLE 11. INTENSITY AND POLARIZATION FOR TITANITE (l 
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TABLE 9. INTENSITY AND POLARIZATION FOR QUARTZ (l 
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TABLE 8. INTENSITY AND POLARIZATION FOR QUARTZ (l 
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INTENSITY AND POLARIZATION FOR QUARTZ (l 
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TABLE 6. INTENSITY AND POLARIZATION FOR QUARTZ (l 
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TABLE 5. INTENSITY AND POLARIZATION 
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TABLE 4. INTENSITY AND POLARIZATION FOR ICE (l 
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TABLE 2. INTENSITY AND POLARIZATION FOR ICE (l 
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TABLE 1. INTENSITY AND POLARIZATION FOR ICE (l 
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TABLES AND GRAPHS 
FOR POLARIZATION AND INTENSITIES 


In Tables 1 through 40 the normalized intensity, the total normalized 
intensity, and the degree of polarization are given for each 5 degrees of the 
angle for the following materials: 


Water (ice) 

m = 1.33 

- O.OOi 

Water (dirty ice I) 

m = 1. 33 

- 0. 05i 

Water (dirty ice n) 

m = 1. 33 

- O.lOi 

Titanite 

m = 2. 00 

- O.OOi 

Quartz 

m = 1. 55 

- O.OOi 

Iron 

m = 1. 27 

- 1.371 

Zinc 

m = 1.41 

- 4. lOi 

Graphite 

m = 1. 59 

- 0. 66i 


The particle size ranges are l(5)25. 

Only the total intensity and degree of polarization are shown in Figures 
2 through 81. In the tables: 


XII INTEN 
XI2 INTEN 


= ll(0) 

= I 2 (0) 


INTENSITY 


= -(XI1 + XI2) 


POLARIZATION = 


XII - XI2 
XII + XI2 


THETA 


= angle between the incident 
and scattered direction. 


The'heading Gaussian Integral should be ignored since it is an artifact 
of the programming and contributes no information to the tables. 
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For the Gauss-Legendre quadrature, the are the s zeros of the 

Legendre polynomials P (x), and the W are weighting functions. Both W 

ns s 

and y are tabulated [ 6 1 . In the trapezoidul quadrature, h is the increment 
s 

of integration. It has been found that the size of the integration increment 
is extremely important f 7 1 . To study this effect, a 5-, 16-, and 96-point 
Gauss-Legendre scheme and the trapezoidal quadrature, for comparison, 
were used separately to calculate ij, i 2 , and the polarization. A discussion 
of the results is given in the appendix. 


Distribution 

The only distribution function included in this report is the inverse 
power law: 

n(cv) = b« - , 

with k = 4. It is limited to the inverse power in order that n (a) decreases 
as a increases. This distribution function is useful when using a range of 
sizes that has a certain decay rate. The term b is related to the total number 
of particles per unit volume N' by 

, ( k + 1) N* 

b ~ -k + 1 -k + 1 ’ 

a- . -a 
mm max 

which is found simply by integration. The integral diverges at O' = 0, and, 
hence, a finite lower limit must always be included. The slope of the 
In n(o! ) _versus-ln a curve is easily found and is 

d In n(<v) 
d In (a) ~ 

Figure 82 shows the range and values of the distribution for several 
values of k. Also, it should be noted that regardless of the power of k , the 
smaller size particles dominate the inverse power distribution. 

The distribution function in terms of the radius r is given by 
n(r) ^ b (^) r ' k • ; ■ ■ ■ 
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To determine polarization effects, the degree of polarization defined 
by 

p/ fl \ = i,(m, a?, 6) - io(m, o>, 9) 
i](m, a, 6) + i 2 (m, a, 9) 

is used. For the general case of the nonhomogeneous polydispersed cloud, 
the expression for P-(0 ) becomes- • ■ 


P(0) = 


Q a' . ■. . Q oi 

Yj j maX i 1 (m , a, 0.)n(a')da - J maX i 2 (m , a, 0)n(G)d< 

q=l a . q q=l a . q 

_ mm mm 

Q a 

/. f maX [ii(m a, 0)n(«)dG + i ? (m , a, 0)n(a;)]da 

— I , J q- . q 

■ q=l a . 

mm 


LG 


Two methods have been used to program these equations for the 
computer, a Gauss-Legendre quadrature and the trapezoidal quadrature. The 
integrals were normalized by the total number of particles in a unit volume 

a • ’ g 

I ■ max 

■> n(a ) d« 

a 

mm 

and it is this normalized scattering function which is given in Tables 1 through 
40 and Figures 2 through 81. With the Gauss-Legendre scheme [4] 

b S 

J i (m, g, 0)n(G)da' = — Yj W g i (m, 6 , y g ) n (y g ) > 
a s=l 


where 

The trapezoidal quadrature is 


b 

J f(x)dx 
a 


^ ♦ f(a + 


1-) 


+ f (b - 1) + 


£(by 

2 
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and 


I (0) 

total 



If the cloud is monodispersed and homogeneous, then 


I (0) 

total 



a, 0 ) + j., (m, at, 0) 
2 



where N is the total number of particles in the scattering volume. For a 
homogeneous polydispersed collection of particles 


L i t (m, a , 0 ) + i 2 (m, a , 0 ) 

V n(a ) Aa , 

£=1 * 1 

where n(o! ) is the distribution function for the particle sizes. The term 
n(a ) plays the part of a weighting function. The total number of particles 

per unit volume is given by ^ n(«) A a ; hence N = x volume. 

If the distribution is smooth and continuous, then the summation can be 
replaced by an integral 


I (0) = 
total 


I (01 

total 


and 


2 


a 

r max r . , . 

J [ij (m, a , 0) 

Oi . 

min 


+ % (m, a , 0)] n(a) da 


a 

/ TXlclX 

i (m, a, 0)n(a) da , where c = 1 , 2 . 
c 

a . 
mm 

In addition, if the cloud is nonhomogeneous but consists of several 
different materials, then 


1 (0) \ J maX [ii(m , a, 0 ) + ^ (m , a, 0 )Jn (a)da . 

total q=l & a q q 

min 


4 


A 

e 

j 

n+% 

„< 2 > 

n+V> 

The terms and 77 -S give the radiant intensities of the scattered 
87r 87T & 

components polarized perpendicular and parallel, respectively, to the plane of 
scattering when the sphere is illuminated by unpolarized monochromatic inci- 
dent light of unit intensity. The plane of scattering is defined as the plane 
containing the sphere, the incident light direction, and the scattered light 
direction. (See Fig. 1. ) 

The differential scattering cross section is the ratio of the energy 
scattered in a given direction per unit solid angle per unit time to the energy 
incident per unit area per unit time. It is given by Penndorf [5] . 


= wavelength of the incident light, 

= scattering angle, 

= mev , 

= half-integral order Bessel function, 

= half-integral order Hankel function of second kind. 


A 2 L X 2 i„ 


do - 

dft 


X 2 

87 ? 


(ii + h) 


(area/steradian) 


Some authors call this the angular scattering c ( ross section. 

The Mie theory can easily be extended to a cloud of particles if the 
light reaching the point of observation has encountered only one scatterer and 
if the particles are randomly oriented in space so that there are no phase 
effects. With these assumptions, the total scattered intensity in a particular 
direction is the sum of the intensities scattered from each particle, 


U 0 ) = 2 Z * (“>“(> e ) where c=l, 2 , 

q = 1 1 = 1 4 


3 


and 


i 2 (m, a, Q ) = 


l 

n=l 


dP (x) 
n 

dx 


(1 - x) 2 


d’P (x)' 


n 

dx 2 


+ B 


dP (x) 
n 

n dx 


where 

ij and L, are scattered intensity functions, 
P n (x) is a Legendre polynomial of order n, 
x = -cos 6 . 


and are amplitude functions given by 

dip (P) dip (a) 

, = /.2n+l\ 2n + 1 ' m ^n (g) 

n V 1 / n(n + l) dip (j3 J d£ (a ) 

f(a) — - m ip (/?) — ^ 

d/3 n da 


and 


B = 
n 


d ip (a ) di P n (P) 

/ .2n+l\ 2n + 1 ^n^ da ’ m ^ n ( ff ) <p 

\ 1 / n(n + l) d£ (a) dip ((3) 




where 


U- 


ip ' = Ricatti Bessel function 


ioa(^) J (o) , 

V ’ n+% 


n 


l / 2 

= Ricatti Bessel function (~\ (a) , 

V 2 ' n + *4 


a = 


2ira 


= — — = particle size parameter , 
A 


m = complex refractive index , 


I 
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TECHNICAL MEMORANDUM X 64674 


A CATALOGUE OF NORMALIZED INTENSITY FUNCTIONS 
AND POLARIZATION FROM A CLOUD OF PARTICLES 
WITH A SIZE DISTRIBUTION OF a-* 

INTRODUCTION 

The so-called theory of light scattering from single spheres can be 
extended to a collection of spheres of different sizes. The sizes will vary 
according to some distribution function. Using the distribution function, one 
can calculate the average scattering cross section for any angle. If such cross 
sections are calculated beforehand, the computer time required for light 
scattering calculations can be shortened. It is with this thought in mind that 
this report is written. The materials have been limited primarily to those 
that have been used in modeling the astronomical phenomenon, zodiacal 
light, where single particle scattering can be considered sufficient. 

In the course of studying the effects of distributions, it became clear 
that certain parameters strongly influence the results of the scattering by a 
distribution of particles. They are the obvious physical parameters, particle 
size range, refractive index, and the not-so-obvious numerical parameter, 
the integration increment. The effects of all three parameters are investigated. 

Only particles of spherical shape are considered since this is the only 
shape for which the problem of scattering of electromagnetic radiation has 
been solved exactly in closed form and which is easily adapted to the computer. 


THEORY 

The development of the Mie theory is repeated in the literature [ 1-3] . 
The basic expressions [4] for the radiation scattered from a sphere of radius 
a and refractive index m are: 


i t (m, a, 9 ) = 


CO 

V 

u 

n=l 


dP (x) 
n 

dx 


+ B 


n 


dP (x) 
n 

dx 


- ( 1 - x )* 


d 2 P n (x)| 

dx 2 
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